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Abstract 

Models for supernovae (SNe) arising from thermonuclear explosions of white dwarfs (WDs) have been 
extensively studied over the last few decades, mostly focusing on the single degenerate (accretion of 
material of a WD) and double degenerate (WD-WD merger) scenarios. In recent years it was suggested 
that WD-WD direct collisions provide an additional channel for such explosions. Here we extend 
the studies of such explosions, and explore the role of Helium-shells in affecting the thermonuclear 
explosions. We study both the impact of low-mass helium 0.01 M 0 ) shells, as well as high mass 
shells (> 0.1 M 0 ). We find that detonation of the massive helium layers precede the detonation 
of the WD Carbon-Oxygen (CO) bulk during the collision and can change the explosive evolution 
and outcomes for the cases of high mass He-shells. In particular, the He-shell detonation propagates 
on the WD surface and inefficiently burns material prior to the CO detonation that later follows 
in the central parts of the WD. Such evolution leads to larger production of intermediate elements, 
producing larger yields of 44 Ti and 48 Cr relative to the pure CO-CO WD collisions. Collisions of WDs 
with a low-mass He-shell do not give rise to helium detonation, but helium burning does precede the 
CO bulk detonation. Such collisions produce a high velocity, low-mass of ejected burned material 
enriched with intermediate elements, with smaller changes to the overall explosion outcomes. The 
various effects arising from the contribution of low/high mass He layers change the kinematics and 
the morphological structure of collision-induced SNe and may thereby provide unique observational 
signatures for such SNe, and play a role in the chemical enrichment of galaxies and the production of 
intermediate elements and positrons from their longer-term decay. 


1. INTRODUCTION 

Thermonuclear explosions of white dwarfs (WDs) have 
been extensively studied over the last decades, and are 
thought to be the progenitors of various types of super¬ 
novae (SNe), most notably type la SNe, but p ossibl y 
also calcium-rich faint type lb SNe (jPerets et al.l 12010 ). 
Models for thermonuclear detonations are traditionally 
divided between double degenerate (DD) models and sin¬ 
gle degenerate models (SD). The former focus on cases of 
mergers of two WDs, and the latter on accretion o£ma- 
terial from a companion star onto a WD (See lMaoz et al.1 
120141 for a review). In recent years it was suggested that 
the DD scenario can be extended to WD-WD collisions 
and not only mergers, but the former were thought to 
be extremely rare, and occur only in dense stellar clus¬ 
ters. For this reason they attracted relatively little at¬ 
tention compared with other WD explosion progenitors. 
Such collisions , however, are likely to be observable as 
type la SNe (Rosswog et al. 200*1 ) and possibly non¬ 
standard SNe (jRaskin et al.l 120091 . 120101 ). Recently, it 
was shown that some triple systems may dynamically 
evolve through a quasi-s ecular process, reminisce nt of 
Kozai-Lidov oscillations (jKozail 1 19621 : iLidovI [T962) , but 
where significant p eri-center changes can occ ur on a sin¬ 
gle orbit time-scale lAntonini fe Peretsl (j2012D . leading to 
extremely close peri-center approaches. In particular, it 
was suggested that such evolution in triples hosting an 
inner WD-WD binary could lead to ph ysica l coll isions 
and the production of type la SNe (jTdLompsonl 120111: 
iKatz fe Pond [20121 : IKushnir et al.ll2013D . Though such 
system s are not likely t o be frequent (e.g. iHamers et al.l 
( 2013 ): see also section l4~3l in the following), they may 


be much more frequent than WD-WD collisions in dense 
clusters, making such scenarios more likely. 

While previous studies mostly focused on collisions of 
carbon-oxygen (CO), the potential role of a He compo¬ 
nent, and in particul ar He-shells on CO-WD s attracted 
only little attention. IKushnir fe Katzl ( 2014 ) calculated 
the conditions for a detonation in a tiny helium shell 
around two 1M 0 WDs and concluded that detonation 
cannot occur during the collision. Realistic CO WDs typ¬ 
ically hold outer Helium layers which can range in mass 
between 10 -3 M 0 for more mas sive WDs up to 0.024 
M^ f or the lowest mass CO WDs (lLawlor fe MacDonaldl 
l2006f h In the case of an interacting binary where a 
WD grows through accretion, the accreting WD might 
have a significant outer Helium l ayer of up to 0.1 M 0 
or more on low mass WD (e.g. Bildste n et al.1 (2007): 
such a helium layer is thought to play a major role in the 
SD sub-Ch andrasekhar and .la models for thermonuclear 
explosions jW ooslev et al. 1986; Livne fe G lasner 19911: 

Bilds ten et al. 1200 7: Shen & Bildsten 2009: IPerets et all 

201d : IWaldman et al.ll201 lh : the helium layer may ignite 

and later possibly induce the deto nation of the CO-WD 
bulk through converging shock (jShen fe Bi ldsten 2014). 
Recently, it was suggested t o play a similar role in WD- 
WD mergers iPakmor et al.l ( 201 3). Here we study the 
outcomes of WD-WD collisions and explore the possibil¬ 
ity that a helium shell can play an important role in such 
collisions. As we discuss in the following, a He-shell can 
produce non-negligible effects on the collisions outcome, 
changing the composition, energetics and velocity distri¬ 
butions of the resulting SN; the amplitude of these effects 
depend on the mass of the He layer considered, with more 
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significant effects arising from larger He layers. 

We begin by describing our 2D simulations and their 
initial conditions (section^), and present the main results 
in section [3l We then discuss our results (section 0]), 
where we also discuss the rates of such collisions and 
point out their main implications and their observational 
signatures. 

2. SIMULATIONS 

2.1. WD Models 

The WD models are calculated using the MESA stel¬ 
lar evolution code dPaxton eDaHl2Qlll ). The initial WD 
models were taken from IWolf et ahi (|2013[ ) and then 
changed using MESA for the different compositions. The 
models with a small helium shell were built by relaxing 
the chemical composition of the non helium elements to 
be 50% carbon and 50% oxygen. For the CO models we 
removed the helium shell up to a total of 10 _3 M©. After 
that we relaxed the WDs composition to be 50% carbon 
and 50% oxygen. The models with a large helium shell 
were built by accretion at a rate of 10 -8 M©yr -1 of he¬ 
lium to the initial models up to a total 0.1 M© — 0.2M© 
of helium on the WDs. In the following we consider the 
whole range of low-mass to high mass helium shells, but 
we should note that the latter cases of massive shells 
are likely to be rare. Nevertheless, even massive helium 
shells may exist around WDs in binaries with a helium 
donor star (e.g. AM-CVn systems), and the possibility of 
helium shell detonation in such cases has been widely ex¬ 
plored (e.g. IWooslev et ahl 1 19861 : iLivne fe Glasnerl 1 19911 : 

IWaldman et al .112 01 lL and references therein). A collision 
of a CO WD with such a WD with a massive helium shell 
could therefore occur either randomly in dense clusters, 
or in high multiplicity systems (>2 component) in which 
one of the components is such a helium accreting WD 
in a close binary, and another outer CO WD companion 
exists (e.g. llKrzeminskill 19721 find a triple system with an 
inner AM CVn). The outer WD do not need to be close, 
as even collis ions with very far comp anions could occur 
in the field (jMichaelv fe Peretsl 120151 ). 

2.2. Hydrodynamical Simulations 

We carried 17 different 2D simulations of WD- 
WD head-on collisions using the FLASH v4.2 code 
(jFrvxell et alJ 12000 ). The widely used FLASH code is 
a publicly available code for supersonic flow suitable 
for astrophysical applications. The simulations were 
done using the unsplit PPM solver of FLASH in 2D 
axisymmetric cylindrical coordinates on a grid of size 
3.2 x 6.4 [lO 10 cm] using adaptive mesh refinement 

(AMR). To prevent the production of unreal early deto¬ 
nation that may arise from nume rical resolutio n we ap - 
plied a limiter approach following iKushnir et al.l ( 2013 ). 
We made multiple simulations with increased resolution 
until convergence was reached in the nuclear burning. We 
found a resolution of 5 — 10 km to be sufficient for con¬ 
vergence to up to 10%. Gravity was included as a mul- 
tipole expansion of up to multipole Z = 10 using the new 
FLASH multipole solver. The equation of state used in 
the si mulations is the Helmholtz EOS (jTimmes fe Swestvl 
2000). This EOS includes contributions form partial de¬ 
generate electrons and positrons, radiation, and non de¬ 
generate ions. The nuclear network was FLASH’S 19 


Model Primary Secondary 

Mass [M 0 ]Pc [g]Mass [M 0 ] p c [g] 


1A 

0.6 

CO 



3.4 

X 

10 6 

0.6 

CO 


3.4 

X 

10 6 

IB 

0.6 

CO 



3.4 

X 

10 6 

0.6 

CO 

+ 0.01 He 

3.5 

X 

10 6 

1C 

0.6 

CO 

+ 

0.01 He 

3.5 

X 

10 6 

0.6 

CO 

+0.01 He 

3.5 

X 

10 6 

ID 

0.6 

CO 

+ 

0.024 He 

4.0 

X 

10 6 

0.6 

CO 

+0.024 He 

4.0 

X 

10 6 

IE 

0.6 

CO 



3.4 

X 

10 6 

0.6 

CO 

+ 0.1 He 

5.7 

X 

10 6 

IF 

0.6 

CO 

+ 

0.1 He 

5.7 

X 

10 6 

0.6 

CO 

+ 0.1 He 

5.7 

X 

10 6 

1G 

0.6 

CO 

+ 

0.2 He 

1.1 

X 

10 7 

0.6 

CO 

+ 0.2 He 

1.1 

X 

10 7 

2A 

0.7 

CO 



5.9 

X 

10 6 

0.7 

CO 


5.9 

X 

10 6 

2B 

0.7 

CO 



5.9 

X 

10 6 

0.7 

CO 

+ 0.01 He 

6.1 

X 

10 6 

2C 

0.7 

CO 

+ 

0.01 He 

6.1 

X 

10 6 

0.7 

CO 

+0.01 He 

6.1 

X 

10 6 

2D 

0.7 

CO 

+ 

0.024 He 

6.8 

X 

10 6 

0.7 

CO 

+0.024 He 

6.8 

X 

10 6 

2E 

0.7 

CO 



5.9 

X 

10 6 

0.7 

CO 

+ 0.1 He 

9.9 

X 

10 6 

3 

0.6 

CO 



3.4 

X 

10 6 

0.7 

CO 


5.9 

X 

10 6 

4 

0.7 

CO 



5.9 

X 

10 6 

0.8 

CO 


1.0 

X 

10 7 

5A 

0.8 

CO 



1.0 

X 

10 7 

0.6 

CO 


3.4 

X 

10 6 

5B 

0.8 

CO 



1.0 

X 

10 7 

0.6 

CO 

+ 0.01 He 

3.5 

X 

10 6 

5C 

0.8 

CO 



1.0 

X 

10 7 

0.6 

CO 

+ 0.1 He 

5.7 

X 

10 6 


TABLE 1 
Collision models. 

isotopes alpha chain network. This network can cap¬ 
ture welLthe en ergy g enerated during the nuclear burning 
(|Timmes et al.lf2000t ). As initial conditions we take the 
two WDs to be at contact with free fall velocities. The 
various simulations conducted are summarized in table 

m 

3. RESULTS 

3.1. Collisions of pure CO-WDs and comparison with 
previous collision simulations 

We run several pure CO-WDs collisions in order to 
check the consistency of our simulations with previous 
works, and then use them as a benchmark for comparison 
with the novel collision simulations in which WDs with 
He-layers were considered. 

Previous collision simulations include several works 


by different groups (Rosswog et al. 200jl Raskin et al. 

2009 

. 120101: Loren-Aguilar et al.l |2010|: Hawley et al. 

2012 

: Kushnir et al.| 2013|: Garcfa-Senz et ahl 20131). In 


these simulations only pure CO-WD collisio ns were s im¬ 
ulate d (or CO-WD - He-WD collisions Ros s wo g et al.1 
2009). Early studies of WD-WD collisions showed large 
discrepancies in the 56 Ni yields produc ed, likely result¬ 
ing from insufficient resolution (jKushnir et al.l l2013) and 
the use of too-large time steps in couplin g th e nu¬ 
clear burning to the hydrodynamics (jHawlev et alJf2QT2h . 
IKushnir et al . (IMl devised the limiter approach to 
tackle such potential difficultie s; recent high resolutions 
fe.g.. iGa rcfa-Senz et al. 2013) show a convergence to 
IKushnir et al.l (|2013h results, supporting the validity of 
their limiter approach. Our main com pariso n is there¬ 
fore with the results bv IKushnir et ah (201 3). who sim¬ 
ilarly used the FLASH code, and included the limiter 
approach. One difference in our simulations i sthejuse of 
the unsplit hydro solver in FLASH while IKushnir et al.l 
120131 uses the split solver. 

Our pure CO-WDs collisions result in CO detonations 
in both WDs. The total energy Ep and abundances of 
helium, intermediate elements, and nickel are presented 
in Table 13.21 All the results are shown after the nu¬ 
clear reactions in the simulations ended. includes ki- 
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netic, thermal, and gravitational energy. The amount of 
nickel^and total energy we get is similar to the results of 
iKushnir et al.l (|2013h . This is expected as in both works 
the same code is used with the same limiter on the nu¬ 
clear reaction rate. We find our results to be similar up 
to ^ 10% in the pure CO cases. These small differences 
should exist as a result of the slightly different initial 
preparation of the progenitors and the resolutions used 
in the two studies. We therefore confirm previous stud¬ 
ied of such collisions and find them to be consistent with 
the results we obtain 

3.2. Collisions between CO-WDs with a low-mass, 
non-detonating helium layer 

Collisions of WDs with low-mass helium layers (up to 
Mjie < 0.024M©; see Fig. [Tj) do not give rise to a det¬ 
onation in the helium layer. Such collisions are char¬ 
acterized by three phases: a) The helium compression 
phase. A short time after the collision the helium layer 
is compressed to densities of pn e ~ 10 5 gem -3 and tem¬ 
peratures of The ~ 10 9 K. Under these conditions the 
helium stars burning, but no detonation is triggered. 
The compressed helium is ejected in the plane perpen¬ 
dicular to the collision direction with velocities of up to 
20,000kms -1 as it burns to intermediate elements, b) 
The CO detonation phase. This phase is very similar 
to the pure CO-WDs collisions. The compression and 
burning of the helium layer has almost no effect on the 
CO core. The detonation in the CO material is trig¬ 
gere d in a very similar way to pure CO-CO collisions 
([Kushnir et al.ll20l3l ) and can be seen in Fig. [lj c) Outer 
helium layer ejection. The detonation wave in the CO 
material reaches the outer part of the helium layer and 
eject the partially burned helium layer. At the end of the 
explosion ^ 1 /2 of the initial helium mass survives in all 
cases and is ejected after the shock reaches the outer part 
of the WD. Intermediate mass elements are also ejected 
in this way. The 56 Ni yields in all cases are compara¬ 
ble to the 56 Ni yield in pure CO-CO collisions. The total 
mass of the intermediate elements is summarized in table 
13.21 We caution that exact outcomes of these simulations 
depend on the nucleosynthetic network used. Here we 
use the 19 isotopes alpha chain network; in section 13.41 
we show that although the simplified smaller network 
well captures the overall evolution, energetics and nucle¬ 
osynthetic production, the result obtained differ quanti¬ 
tatively (and non-negligibly) from those obtained using 
a much larger network. In particular, the latter show 
an increase of the nickel production and suppression of 
the IME production, compared with the results of the 
smaller network. 

3.3. Induced helium detonation in collisions between 
CO-WD with a large-mass helium layer 

In collisions of WDs with a massive (Mn e > 0.024M©) 
helium layers, the helium detonates during its compres¬ 
sion phase. For the Mn e = 0.024M© case the helium 
detonation is not able to propagate through the curved 
helium shell (Fig. [2]) and has almost no influence on 
the overall collision outcomes, i.e. the detonation in the 
CO core is triggered by the collision in a similar way 
to pure CO-WDs collisions. For more massive shells 
(Mn e = O.1M 0 ; models 1F,2E) the helium detonation is 


Helium 


Model 

4 He 

40 Ca 

44 T j 

00 

Q 

I-S 

56 Ni 

Eexp 

Detonation 

1A 

l(-3) 

4(-2) 

9(-5) 

4(-4) 

0.31 

1.3 

- 

IB 

5(-3) 

4(-2) 

1.5(-3) 

l(-3) 

0.27 

1.3 

No 

1C 

l(-2) 

5(-2) 

3(-3) 

2(-3) 

0.29 

1.3 

No 

ID 

3(-2) 

4(-2) 

5(-3) 

3(-3) 

0.38 

1.2 

Yes 

IE 

4(-2) 

4(-2) 

1.3(-2) 

l(-2) 

0.4 

1.4 

Yes 

IF 

1(-1) 

4(-2) 

2(-2) 

2(-2) 

0.8 

1.8 

Yes 

1G 

1(-1) 

6(-2) 

7(-3) 

9(-3) 

0.3 

2.1 

Yes 

2A 

3(-3) 

5(-2) 

l(-4) 

6(-4) 

0.52 

1.5 

- 

2B 

6(-3) 

4(-2) 

8 (-4) 

9(-4) 

0.65 

1.7 

No 

2C 

8(-3) 

5(-2) 

1.5(-3) 

l(-3) 

0.54 

1.6 

No 

2D 

3(-2) 

5(-2) 

3(-3) 

3(-3) 

0.61 

1.6 

Yes 

2E 

5(-2) 

5(-2) 

1.4(-2) 

l(-2) 

0.67 

1.6 

Yes 

3 

2(-3) 

4(-2) 

l(-4) 

5(-4) 

0.42 

1.4 

- 

4 

1 (-3) 

4(-2) 

1.6(-4) 

6(-4) 

0.49 

1.6 

- 

5A 

2(-3) 

3(-2) 

4(-5) 

3(-4) 

0.32 

1.4 

- 

5B 

7(-3) 

4(-2) 

l(-3) 

6(-4) 

0.27 

1.4 

No 

5C 

5(-2) 

6(-2) 

1.5(-2) 

l(-2) 

0.4 

1.7 

Yes 


TABLE 2 

Abundances at the end of the simulations in values of 
M 0 . Energy is in 10 51 erg. 


able to propagate through the shell before the CO is det¬ 
onated (see Fig. m. In this case we find that the helium 
shell burning results in the production of non-negligible 
amounts of intermediate mass elements and nickel. The 
velocity of the helium-burned material is lower than the 
velocity of the material ejected following the CO core 
detonation which gradually outruns the material ejected 
through the earlier He-detonation in the shell. For even 
more massive helium shells ((Mn e = 0.2M 0 ; model 1G) 
the converging shock caused by the helium detonation 
is able to detonate the CO. In this case the induced CO 
detonation is obtained earlier than the corresponding CO 
detonation in the pure CO-CO WD collisions and both 

log p [g cm -3 ] 

1 2 3 4 5 6 7 


2 

Td 1 

o 0 
o, 

N -1 

-2 
2 

o 0 
6 

N -1 

-2 
2 

o 0 
o, 

tsi -1 
-2 

012340123401234 
X [0.01fl Q ] X [O.Olfl 0 ] X [0.0LRJ 

Fig. 1. — Density before,during, and after CO detonation for 
two O.6M 0 WDs. Upper panels: pure CO WDs. Middle panels: 
the top WD with O.OIMq of helium while the bottom WD is a 
pure CO. Bottom panels: both WDs with a layer of O.OIMq of 
helium. The detonation in the CO is similar in the three cases. 
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Fig. 2.— Ignition in the He layer and C/O core of two WDs 
with composition of 0.6Mq C/O and 0.024 Mq of He. The helium 
layer is detonated first resulting in a fast ejecta in the horizontal 
plan. The helium detonation front is not able to propagate along 
the curved helium shell. The helium ignition has no real effect on 
the CO ignition process. 


lQ g PTm 

-6 -4 -2 0 2 4 



012345601234567 
X [0.01i? Q ] X [0.01i? Q ] 


Fig. 3. — Density of 44 Ti after the shock reaches the outer parts 
of the WDs. A) Top WD is a O. 6 M 0 CO with 0.01 Mq of helium; 
bottom WD is a pure 0 . 6 Mq CO WD. B) Top WD is a 0.7MqCO 
with O.OIMq of helium; bottom WD is a pure 0.7Mq CO WD . 
C) Both WDs are 0.6MqCO with O.OIMq of helium. D) Both 
WDs are 0.6MqCO with 0.024Mq of helium. Most of the ejected 
44 Ti is from the outer parts of the WDs from partial burning of 
the helium layer. A small plump is ejected during the collision and 
can be seen in front of the WDs. 


the nickel abundance (M Ni se ~ 0.3M 0 ) and the abun¬ 
dance of intermediate elements are smaller than in the 
pure CO collisions. 

3.4. Nuclear reactions post-processing 

In order to provide a better understanding of the nucle- 
osynthetic yields from WD-WD collisions we go beyond 
previous models using a small nuclear reaction netwe- 
orks networks. We extend these by including a large 160 
post-processed nuclear reaction network for two of our 
simulations, a pure CO-CO collision ( 1 A) and a CO+He - 
CO+He (low-mass helium shell) collision (1C). The post¬ 
processing was done by including 30,000 trace particles to 
the simulations and record their thermodynamical con¬ 
ditions throughout the simulations. We post-processed 
the results from the trace particles using the burn unit 
of the MESA code (version 7624) ([Paxton et al.l 2015 ). 
MESA uses the J IN A reaclib version V2.0 2013-04-02 
(iCvburt et al.llll 2010 ) database for nuclear reaction rates. 
The results are summarized in table [3j 

The results show production of a larger amount of 
56 Ni and suppression of 44 Ti with respect to the re¬ 
sults obtained using the smaller a network in FLASH 
(19 . Similar differences due to the use of larger nu- 
cleosynthetic network h ave been seen in other works 
(jGarcfa-Senz et all 120131 ). which arise due to the more 
efficient nuclear burning, producing larger abundances 
of higher elements (e.g., 56 Ni) and lower amounts of in¬ 
termediate elements (e.g., 44 Ti). 

3.5. Influence of initial separation of the WDs 

The simulations in this study were initiated when the 
two WDs are at contact, and the effect of tidal distor¬ 
tion of the WDs prior to contact were not taken into 
account. To verify whether such pre-impact effect play a 
role in the collision outcomes we ran an additional simu¬ 
lation, identical to 1 C, but with an initial separation of 
3 (.Rwdi + i^wD 2 )- The results show vary little compares 
with those obtained from ‘initial contact ‘ simulations, in 
particular there only a very small distortion is observed 
in the WDs structure prior to impact. This is antici¬ 
pated as the free fall velocity is supersonic and the WDs 
do not have sufficient time to be significantly distorted 
by much before the collision. In summary, the energetics, 
overall evolution post impact and nucleosynthetic yields 
show only very slight deviations from the ‘initial contact‘ 
simulations. 

4. DISCUSSION 

4.1. Composition, morphology and kinematics 

In our simulations of WD-WD collisions including a 
He-layer we find that the helium is burned mostly into 
intermediate elements resulting in a much higher mass 
ejecta of these elements compared with pure CO colli¬ 
sions. Only a small fraction of the helium shell burns 
into 56 Ni, which has only a minor effect on the 56 Ni 
yields. The deviation in the yield of intermediate el¬ 
ements compared with their yield in pure CO colli¬ 
sions is found to be almost linearly dependent on the 
He mass (Fig. [5]). Using LS fit we find the rela¬ 
tion Mhe 4 = 0.001883 + 0.462991M s heii, he., almost 
half of the helium is left unburned. For 44 Ti and 48 Cr 

we find M t *44 = (0.132099 ± 0.006176) M s h e ii, M cr48 = 
(0.094573 db 0.003914) M sheU . 
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0.00.40.81.2 0.00.40.81.2 0.00.40.81.2 0.00.40.81.2 
X [0.01i? o ] X [O.Oli? 0 ] X [0.017? Q ] X [O.Oli? 0 ] 


Fig. 4.— Ignition in the Helium shell and the CO core of two WDs with composition of 0.6Mq CO and 0.1 Mq of He. The helium layer 
detonates first and propagate along the outer shell. The CO core is then ignited as a result of the collision’s shock. 


0.016 


0.014 


0.012 
J 0.010 
j= 0.008 

■o 

cu 

& 0.006 
nr 

0.004 

0.002 


0.00 0.02 0.04 0.06 0.08 0.10 0.12 

Hellium shell mass [M e ] 

Fig. 5. — Mass ejecta of 44 Ti and 48 Cr for the different models. 
The lower left points are for pure CO-WD collisions. 

For helium shells with mass M s h e ii < 1O _3 M 0 these 
yields are significantly higher than those produced in col¬ 
lisions of bare CO white dwarfs. For low mass shells the 
helium ejecta mostly burns as a result of the CO deto¬ 
nation shock and the ejecta is relatively spherically sym¬ 
metric besides a small ejected plume of He/He-burned 
material (see Fig. n. However, a bipolar structure is 
formed when the two helium shells on the corresponding 
two WDs have different masses (Fig [3]). In the high mass 
cases most of the burning results from the helium deto¬ 
nation wave. A large fraction of intermediate elements 
are then ejected in a cone of « 30° to the collision plane. 

The velocity distribution of the various components 
also differs between the different type of collisions. In 
Fig. [6] we present the velocity profiles for cases 1A, IB, 
IF in the upper, middle, and lower parts of the simu¬ 
lation domain. Collisions which include He-layers show 
high velocity tails of ejecta of intermediate elements. In 
particular the most He-rich collisions give rise to a non- 
negligible mass of material ejected at higher velocities 
due to the He-detonation stage. Using a much larger nu¬ 
clear network (see section [3~4|) we find a reduction of IME 
production and an increase of nickel production, show¬ 
ing that a large nuclear network is re quired for p roduc¬ 
ing the correct abundances (see also lHolcomb fc Kushnirl 
120151 for pointing out these issues in a similar context). 



4.2. Comparison with other models of type la 

The overall energetics and behavior of the WD-collision 
models can be relatively similar to other models for type 
la SNe. However, high yields of intermediate elements 
do not show in models of t he single degen e rate s cenario. 
For example, model W7 of iNomoto et al.l (j 19841 ) results 
with < lO _4 M 0 of 44 Ti and < 1O -5 M 0 of 48 Cr. As the 
amount of helium on a Chandrasekhar mass WD is small 
we do not expect high yields of these eleme nt^ in an y case 
except maybe for a failed l a supernova (e.g.. I Jordanet al.l 
120121 : iKromer et al.l [2013: this resu lt might also is ap¬ 
ply to the core-dege nerate scenario iLivio fc Riessl 120031 
iKashi fc Sokerll201l[ ). 

In the double deg enerate scena rio a helium shell can 
exist on the WDs. [Raskin et al.l 1 2012 ) simulated the 
merger of two WDs with a helium layer; they found 
that a detonation occurs in the helium layer when the 
primary mass was 1.06M 0 ; however it did not lead to 
car bon deto nation as a result of the helium detonation. 
iPakmor et al. ( 201 3) found that a violent merger with a 
O.OlM 0 helium shell produces less than 2 x 1O _8 M 0 of 
elements more massive than calcium in the helium shell. 

The relatively large abundances of intermediate el¬ 
ements such as 44 Ti can provide a more unique 
signature for He-rich collisions, though high abun¬ 
dances of these components can also be produced 
in the SD sub-Chandrasek h ar and .la explosion 
models (IWooslev et al.l 119861: iLivne fc Glasnerl f l991: 

Bildsten et al.l 1200 7: Shen fc Bildstenl 12009: Pe rets et al l 

2010; Waldman et al. 2011)). Interestingly, recent obser¬ 
vations (iTroia et al.l 2(3144 suggest that the 44 Ti abun¬ 
dances in Tycho SN remnant are comparable with the 
abundances we find in the Helium rich colli s ions ( even 
those with low-mass He-layer). iLopez et al.l ( 2015 ) did 
not detect any 44 Ti but only derived upper limits of 
2.4 x 10 -4 M 0 . Such 44 Ti levels are significantly higher 
than expected in typical SD scenarios (and much larger 
than the W7 model), but are still consistent with the 
main models (1B,2B); and are lower than the sub- 
Chandrasekhar models (depending on th e as sumed dis- 
tance of the Tycho SNR; see figure 4 in lTroia et al.ll2014l ). 
If 44 Ti rich SNe are frequent, then their large production 
of 44 Ti could also be important for the 511 keV emission 
in the Galaxy due to positron-electron annihilation fol- 
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1000 km s -1 


Fig. 6. — Velocity profiles of 40 Ca, 44 Ti, 48 Cr, 56 Ni from three representing simulations. The simulations are for the cases: 1A - no 
helium, IB - one WD with 0.01 Mq of helium, IF - two WDs with O.IMq of helium. Left: lower region with cos 9 < —1/3. Middle: middle 
region with 1/3 < cos 6 < 1/3 . Right: upper region with cos 6 >1/3. The angle 0 is the polar angle from the symmetry axis. 


lowing 44 Ti decay (see lPeretslll2014j for additional discus¬ 
sion of He-rich SNe contribution). 

4.3. The rate of WD-WD collisions 

The rate of WD-WD collisions in dense stellar systems 
is low, and the overall rate was suggested to be dom¬ 
inated by possible WD-WD collisions in quasi-secular 
triples. Estimating the rate of WD-WD collisions is 
complicated and involves many uncertainties, both in 
the initial fractions and orbital parameters of triple sys¬ 
tems as well as the effects of dynamical and stellar evo¬ 
lution. Such studies are beyond the scope of this pa¬ 
per which focuses on the properties of the explosions 
themselves. Nevertheless, one can constrain an upper 
l imit o n such rat es through simple arguments (see also 
llSoker et 151120141 ). 

Let us consider the most optimistic case. If we assume 
the binary components in wide binaries are formed inde¬ 
pendently (uncorrelated) then the mass-function of each 
stellar companion is independent. In such a case, the 
fraction of WD companions to WDs among WD-binaries 
is of the order of the overall WD fraction among sin¬ 
gle stars, namely ^ 10%. If all of these were members 
of quasi-secular triple systems, with isotropic distribu¬ 
tion of mutual inclinations between the inner and outer 
binaries of the triples, then ^ 5% of these inner WD- 
WD bina ries will co llide due to quasi-secular evolution 
(|Katz fc Pond [2Q12h . with about half of the collisions 
having impact parameter sufficiently small (smaller than 
half the WD radius), as to lead to a detonation and a full 


thermonuclear explosion. In other words, 0.25% of WD 
systems will produce a thermonuclear explosion. Given 
that ~ 1% of WD are required to explain the rates of 
type la SNe, then in the most optimistic and highly sim¬ 
plistic case WD collisions in triples can explain 25 % of 
type la SNe. 

In practice, current available data suggest the rates are 
significantly smaller. Among 0.9 — 1.5 M 0 stars (consti¬ 
tuting half of the potential WD progenitors in the mass 
range 0.9 — 8 M 0 ), the triple and higher multiplicity frac - 
tion is only 13% (|Tokovininl 12014) : [Riddle et al.l (20151 ). 
Moreover, among these, only a small fraction(< 10% ) 
of the systems have period ratios residing in the quasi¬ 
secular regime ( Pont/Pin < 20). Finally, there is rea¬ 

son to believe that the mutual inclinations of at least 
the sho rter per iod triples (< 100 yrs), are relatively 
small (|Fekellll98lT ), rather than having an isotropic dis¬ 
tribution. Such low inclination triples will not con¬ 
tribute at all to the quasi-secular regime. Taken to¬ 
gether, one may conclude that WD progenitor systems 
in lower mass range may give rise to ^ 100 times smaller 
rates than the most optimistic estimated rate considered 
above. Even allowing for a few times higher fractions 
due to some unaccounted-for incompleteness, these ob¬ 
servational data already suggest very low rates. More¬ 
over, some of the potential triple progenitors will induce 
a merger/collision of the inner binary before its evolution 
into a WD-WD binary. More massi ve stars are known 
to have higher multip licity fractions (Duchene fc Krausl 
120131 : [Sana et al.ll2013D . and the rates might therefore be 
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higher in the upper half of the mass range. However, 
such transition in the multiplicity and orbital properties 
with mass is likely to be continuous, and not produce or¬ 
ders of magnitude higher fractions of quasi-secular triples 
than those seen for the lower mass triples. 

We conclude that the rate of WD-WD collisions is 
likely to give rise to only a small fraction of the type 
la SNe rate, of the order of at most a few percents of the 
total type la SN rate in the more optimistic case, and 
likely lower. Note, however, that the current data on 
triple and higher multiplicity systems are still very lim¬ 
ited, and more data are needed in order provide better 
estimates, especially for collisions of more massive WDs. 

5. SUMMARY 

We investigated the result of head-on collisions of CO 
white dwarfs with a helium layers as a channel for type 
la supernovae. Pure CO WD-WD collisions were sug¬ 
gested in the pa st to be a channel for for type la in 
globular clusters (iBenz et al.lll989h iRosswog et al.ll2QQ9h 
and more recently suggested as a channel for_the_j)ro- 
duction of type la SNe fr om triple systems (jThompsonl 
1201 ll : UTatz fc Dongll2012f h Simulations of pure CO-WD 
collisions could have observational properties compati¬ 
ble with observed type la SNe in terms of the obser ved 
Iron elements abundances and energetics (Kus hnir et al.1 
l20l3l : ICarcia-Senz et aDlMsh . 

In this work we extend these studies to include WDs 
with a helium layer on at least one of the WDs. A low 
mass helium shell of mass up to Mn e < O.O24M 0 should 
exist in WDs according to current stellar evolution mod¬ 
els. More massive helium shells can occur as a result 
of mass accretion from a companion. In these cases a 
quadruple system is probably needed to lead to a col¬ 
lision (as He-rich layer due to accretion requires a very 
close stellar companion). 

We studied WD-WD collisions using the FLASH code, 
exploring a range of masses for both the WD CO bulk 
and for the outer Helium layers. We find that collisions 
involving low mass He-layers Mn e = O.O1M 0 give rise 
to He burning but no detonation occurs in the helium 
shell. The helium shell partially burns into intermediate 
elements resulting in a fast ejecta in the perpendicular 
plane to the collision. However, the He-burning does not 
affect the overall evolution of the collision and the CO- 
core detonation due to shock compression. In particular, 
the nickel yields are very similar to those obtained from 
pure CO-WD collisions. Nevertheless, high velocity ma¬ 
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terial rich with 44 Ti, and 48 Cr are produced at an order 
of magnitude higher levels than in pure CO-WD collision. 
We note that in general even pure CO WD-WD collisions 
produce more of such intermediate elements than single¬ 
degenerate models such as W7; the He-enriched collisions 
therefore produce up to two orders of magnitude higher 
levels of such intermediate elements compared with the 
single degenerate case. 

Collisions with more massive He-layers (Mn e = 
O.O24M 0 ) due give rise to He-detonation, but the detona¬ 
tion front does not propagate much before the CO-core 
detonates due to shock compression, and the eventual 
outcomes of the collisions are still very similar to those 
obtained for of Mn e = O.O1M 0 cases. 

When the more He-rich collisions are considered 
(Mn e = 0.1M e ; models 1F,2E) the helium layer is suf¬ 
ficiently thick as to to enable the helium detonation 
wave to propagate along the WD helium-shell. The he¬ 
lium shell is burned mostly to intermediate elements and 
ejected with lower velocities before a detonation is trig¬ 
gered in the CO core. The CO detonation itself is still 
triggered by the shock compression and behaves similarly 
to the case of pure CO-WD collisions, with almost no ef¬ 
fect of the helium burning itself on the CO detonation. 
The nickel yields in this case are also similar to that of 
pure CO-WDs collisions. The CO explosion ejects ma¬ 
terial with velocities higher than that in ejected burned 
helium shell. This results with a second shock wave ac¬ 
celerating the burned helium material to higher veloci¬ 
ties. Finally, in the most He-rich cases (Mn e = O.2M 0 ; 
model 1G), a converging shock from the helium deto¬ 
nation is able to detonate the CO-core earlier, before 
the shock compression due to the collision, resulting in a 
lower production of Iron elements. 

We conclude that WD-WD collisions with He-layers 
may affect the outcome of WD-WD collisions and give 
rise to significant composition changes, as well as possi¬ 
ble kinematic and morphological changes, which can be¬ 
come significant for the most He-rich WD-WD collisions. 
We n o te that after the su bmission of the current paper 
[Holcomb fc Kushnirl ( 2015 ) studied these issues, finding 
similar outcomes and further confirming our results. 
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lABLb 3 Post-processing nucleosyntnetic results tor models 1A and 1C. 

The first column for each model is at the end of the simulation. The second 
column for each model are the abundances after 10 Gy. All results are in 

solar mass. 
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